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Artifact correction 



5 The present invention relates to the field of image processing, for 

example in medical applications. In particular, the present invention relates to a method 
of artifact correction in a data set of an object of interest, to data processing devices and 
to respective computer programs. 

In a CT scanner with a polychromatic source of radiation, such as a 

10 polychromatic x-ray source, a polychromatic x-ray beam passes through matter and 
low-energy photons are absorbed as the linear attenuation coefficient generally 
decreases with energy. As a result, the beam gradually becomes harder, i.e. its mean 
energy increases. The harder the beam, the less it is further attenuated. Therefore, the 
total attenuation is no longer a linear function of the thickness of the absorber. This 

15 non-linear behavior generally leads to a reduction in the attenuation coefficient. In bone 
for example, beam-hardening can cause reductions of up to 10 %. Neglecting this effect 
in the reconstruction process leads to well-known beam-hardening artifacts, such as 
cupping, streaks and flairs. 

The reconstruction of CT images can be done with iterative algorithms. 

20 These algorithms, in particular iterative maximum-likelihood (ML) reconstruction 

algorithms, introduce artifacts along lines of high gradients of absorption values. This is 
well-known from the expectation maximation maximum-likelihood algorithm (EM-ML 
algorithm). Theses artifacts however are not desired, since they lead to a reduction of 
image quality. 

25 

It is an object of the present invention to provide for an improved image 

quality. 

According to an exemplary embodiment of the present invention as set 
30 forth in claim 1, the above object may be solved by a method of artifact correction in a 
data set of an object of interest, wherein an image of the object of interest is 
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reconstructed on the basis of the data set and wherein a statistical weighing is 
performed during reconstruction of the image. 

Advantageously, this may allow for a correction of artifacts in an image 
on the basis of a statistical weighing during reconstruction of the image. By weighting 
5 weighing certain data in the data set statistically the image quality may be improved 
significantly. 

According to another exemplary embodiment of the present invention as 
set forth in claim 2, the data set is a projection data set acquired by means of a source of 
electromagnetic radiation generating a beam and by means of a radiation detector 

1 0 detecting the beam. 

Advantageously, this may allow for an artifact correction of CT images 
acquired in a CT system, which may be used, for example, in the field of medical 
imaging, and therefore for an improved image quality. 

Another exemplary embodiment of the present invention is set forth in 

15 claim 3, wherein the source of electromagnetic radiation is a polychromatic x-ray 
source. Furthermore, according to an aspect of this exemplary embodiment of the 
present invention, the source moves along a helical path around the object of interest 
and the beam has one of a cone beam geometry and a fan beam geometry. 

The application of polychromatic x-rays may be advantageous, since 

20 polychromatic x-rays are easy to generate and provide for a good image resolution. 

Furthermore, since the geometry of the CT scanner system may be of different designs, 
such as, for example, cone beam or fan beam geometry, a method of an exemplary 
embodiment of the present invention may be applied to a plurality of different scanner 
systems and may not be limited to pure CT scanner systems, but also applicable to 

25 positron emission tomography (PET) or single photon emission computed tomography 
(SPECT). 

According to another exemplary embodiment of the present invention as 
set forth in claim 4, the reconstruction of the image is performed on the basis of an 
iterative algorithm comprising a plurality of update steps until an end criterion has been 
30 fulfilled. 

Advantageously, by performing the reconstruction of the image 
iteratively, mathematical reconstruction methods may be used which comprise 
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transcendental equations which can not be solved exactly. Furthermore, an iterative 
approach may provide for improved image quality, since the iterative loop may be 
repeated a plurality of times until the end criterion is met. 

According to another exemplary embodiment of the present invention as 
5 set forth in claim 5, the iterative algorithm is a maximum-likelihood algorithm and the 
reconstructed image has the highest likelihood. Furthermore, the weighing is performed 
in each update step of the plurality of update steps. 

Advantageously, a maximum-likelihood algorithm allows for 
determination of an equation which defines an image with the highest likelihood and a 
10 weighing in each update step of the plurality of update steps allows for an efficient 
compensation or correction of artifacts. 

Another exemplary embodiment of the present invention is set forth in 
claim 6, wherein a number of detected photons during acquisition of the data set is 
determined and the weighing is based on a statistical error of the number of detected 
15 photons. 

Thus, the statistical error of the photon counts is estimated and taken into 
account during the image reconstruction, which may result in an improved image 
quality with less artifacts. 

According to another exemplary embodiment of the present invention as 
20 set forth in claim 7, an update of an attenuation parameter |Lij n+1 is calculated from the 
attenuation parameter |ij n by 

«+l = » . n ^UL ii 

wherein dj is a number of photons emitted by the source of radiation, ly is 
a basis function of an i-th projection, 1* is a vector of basis functions ly of the i-th 
25 projection and </.,//> is an inner product, defined by < l t , ju >= l^Mj • 

This iterative statistical uncertainty weighted algorithm for maximum- 
likelihood transmission tomography may allow for an effective removal of artifacts by 
weighing the update with the intrinsic statistical error of the measured photon counts. 

According to another exemplary embodiment of the present invention as 
30 set forth in claim 8, the reconstruction of the image is based on a sub-set of at least two 
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projections of all acquired projections of the projection data set. 

Advantageously, according to this exemplary embodiment of the present 
invention, performing the statistical weighing on sub-sets of projections may lead to a 
fast convergence of the image with less artifacts. 
5 According to another exemplary embodiment of the present invention as 

set forth in claim 9, a data processing device is provided, which comprises a memory 
for storing a data set and a data processor for performing artifact correction in the data 
set of the object of interest, wherein the data processor is adapted for performing the 
following operation: loading the data set and reconstructing an image of the object of 
10 interest on the basis of the data set, wherein a weighing is performed during 
reconstruction of the image and wherein the weighing is based on statistical 
considerations. 

Advantageously, this may allow for an improved image quality of 
images processed by the data processor which comprise artifacts caused by, for 
15 example, high attenuating objects, such as metal. 

According to another exemplary embodiment of the present invention as 
set forth in claim 10, the reconstruction of the image is performed on the basis of an 
iterative algorithm comprising a plurality of update steps until an end criterion has been 
fulfilled, wherein the iterative algorithm is a maximum-likelihood algorithm and the 
20 reconstructed image has the highest likelihood. Furthermore, the weighing is performed 
in each update step of the plurality of update steps. 

Advantageously, the iterative approach may allow for an improved and 
fast convergence of the image with less artifacts. 

According to another exemplary embodiment of the present invention as 
25 set forth in claim 1 1, a CT scanner system is provided, comprising a memory for storing 
a data set of an object of interest and a data processor for performing artifact correction 
in the data set of the object of interest, wherein the data processor is adapted for 
performing the following operation: loading the data set and reconstructing an image of 
the object of interest on the basis of the data set, wherein a statistical weighing is 
30 performed during reconstruction of the image. 

Advantageously, this may allow for an improved CT image quality of for 
example a part of the body of a patient, wherein the acquired data is processed by the 
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data processor which comprise artifacts caused by, for example, high attenuating 

objects, such as a prosthesis or other implants. 

The present invention relates also to a computer program which may, for 

example, be executed on a processor, such as an image processor. Such a computer 
5 program may be part of, for example, a CT seamier system The computer program, 

according to an exemplary embodiment of the present invention, is set forth in claim 

12. The computer program may be preferably loaded into working memories of data 

processors. The data processors are thus equipped to carry out exemplary embodiments 

of the methods of the present invention. The computer program may be written in any 
10 suitable programming language, for example in C++ and may be stored on a computer 

readable medium, such as a CD-ROM. Also, these computer programs may be available 

from a network, such as the Worldwide Web, from which they may be downloaded into 

image processing units or processors, or any suitable computers. 

It may be seen as the gist of an exemplary embodiment of the present 
1 5 invention that artifacts in an image are corrected on the basis of a statistical weighing 

during the reconstruction of the image. This may lead to an improved image quality 

with less artifacts. 

These and other aspects of the present invention will become apparent 

from and elucidated with reference to the embodiments described hereinafter. 
20 Exemplary embodiment of the present invention will be described in the 

following, with reference to the following drawings: 



Fig. 1 shows a simplified schematic representation of an embodiment of 
25 a computed tomography (CT) scanner according to the present invention. 

Fig. 2 shows reconstructed image slices of a head phantom using the 

maximum-likelihood algorithm (left) and the statistical uncertainty 

weighted maximum-likelihood algorithm according to an exemplary 

embodiment of the present invention (right). 
30 Fig. 3 shows a flow-chart of an exemplary embodiment of a method of 

artifact correction according to the present invention. 

Fig. 4 shows an exemplary embodiment of an image processing device 
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according to the present invention, for executing an exemplary 
embodiment of a method in accordance with the present invention. 



5 Fig. 1 shows an exemplary embodiment of a CT (computed tomography) 

scanner system according to the present invention. With reference to this exemplary 
embodiment, the present invention will be described for the application in medical 
imaging. However, it should be noted that the present invention is not limited to the 
application in the field of medical imaging, but may be used in other applications such 

10 as baggage inspection to detect hazardous materials, such as explosives, in items of 
baggage or other industrial applications such as material testing. 

The scanner depicted in Fig. 1 is a cone beam CT scanner. The CT 
scanner depicted in Fig. 1 comprises a gantry 1, which is rotatable around a rotational 
axis 2. The gantry is driven by means of a motor 3. Reference numeral 4 designates a 

15 source of radiation such as an x-ray source, which, according to an aspect of the present 
invention, emits a polychromatic radiation. 

Reference numeral 5 designates an aperture system which forms the 
radiation beam emitted from the radiation source to a cone-shaped radiation beam 6. 
The cone beam 6 is directed such that it penetrates an object of interest 7 arranged in 

20 the center of the gantry 1, i.e. in an examination region of the CT scanner, and impinges 
onto the detector 8. As may be taken from Fig. 1, the detector 8 is arranged on the 
gantry 1 opposite to the source of radiation 4, such that the surface of the detector 8 is 
covered by the cone beam 6. The detector 8 depicted in Fig. 1 comprises a plurality of 
detector elements. 

25 During a scan of the object of interest 7, the source of radiation 4, the 

aperture system 5 and the detector 8 are rotated along the gantry 1 in the direction 
indicated by arrow 16. For rotation of the gantry 1 with the source of radiation 4, the 
aperture system 5 and the detector 8, the motor 3 is comiected to a motor control unit 
17, which is connected to a calculation unit 18. 

30 In Fig. 1, the object of interest is disposed on a conveyor belt 19. During 

the scan of the object of interest 7, while the gantry 1 rotates around the item of 
baggage 7, the conveyor belt 19 displaces the object of interest 7 along a direction 
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parallel to the rotational axis 2 of the gantry 1 . By this, the object of interest 7 is 
scanned along a helical scan path. The conveyor belt 19 may also be stopped during the 
scans to thereby measure signal slices. Instead of providing a conveyor belt 19, for 
example in medical applications where the object of interest 7 is a patient, a moveable 
5 table is used. However, it should be noted that in all of the described cases it is also 
possible to perform a circular scan, where there is no displacement in a direction 
parallel to the rotational axis 2, but only the rotation of the gantry around the rotational 
axis 2. 

The detector 8 is connected to calculation unit 18. The calculation unit 
10 18 receives the detection result, i.e. the read-outs from the detector elements of the 

detector 8 and determines a scanning result on the basis of these read-outs. The detector 
elements of the detector 8 may be adapted to measure the attenuation caused to the cone 
beam 6 by the object of interest 7. Furthermore, the calculation unit 1 8 communicates 
with the motor control unit 17 in order to coordinate the movement of the gantry 1 with 
1 5 motor 3 and 20 are with the conveyor belt 1 9. 

The calculation unit 18 may be adapted for reconstructing an image from 
read-outs of the detector 8. The image generated by the calculation unit 18 may be 
output to a display (not shown in Fig. 1) via an interface 22. 

The calculation unit which may be realized by a data processor may also 
20 be adapted to perform an artifact correction in the image based on read-outs from the 
detector elements of the detector 8. According to an aspect of the present invention, this 
correction may be performed by reconstructing an image of the object of interest on the 
basis of the read-outs, wherein a statistical weighing is performed during reconstruction 
of the image. 

25 Furthermore, the calculation unit for the data processing device may be 

adapted for performing artifact correction in the data set by performing the 
reconstruction of the image on the basis of an iterative algorithm comprising a plurality 
of update steps until an end criterion has been fulfilled, wherein the iterative algorithm 
is a maximum-likelihood algorithm and wherein the reconstructed image has the 

30 highest likelihood. Furthermore, the weighing is performed in each update step of the 
plurality of update steps. 

Furthermore, as may be taken from Fig. 1, for example, the calculation 
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unit 1 8 may be connected to a loudspeaker 21, for example to automatically output an 
alarm. 

A standard approach for a reconstruction of CT images are iterative 
algorithms, for example a convex maximum likelihood algorithm. The Maximum 
5 Likelihood Method is able to determine an equation which defines an image with the 
highest likelihood L(ju) by setting 8L(jUj)/djUj =0, where ju is the vector of 

attenuation parameters fj,. . Since this equation is transcendental, it cannot be solved 

exactly, but an iterative approach can lead to the solution. One of these iterative 
approaches is the "Convex Algorithm" described in K. Lange and J.A. Fessler, 
10 "Globally convergent algorithms for maximum a posteriori transmission tomography", 
IEEE Trans. Med. Imaging 4, 1430-1450 (1995), which is hereby incorporated by 
reference. The update for one iterative step is found to be 

Y / \de~ <li ^ n> ~y\ 
U n+1 = v" + ju n — ' Vl 1 — , (equation 1) 

J J '?j,<i t .tf 

where d t and Y t are the emitted and observed number of photon counts, 

15 respectively, is the vector of basis functions l i} of the i th projection, and <l n pi> 

is the inner product ^.l {J Mj • This update ensures an increasing likelihood in the 

neighborhood of the optimal point, thus it should lead to an optimal image. 

The reconstruction result of such an iterative algorithm is shown as 
image slice 31 in Fig. 2. Image slice 31 is an image slice of reconstructed projection 
20 data of a forbild head phantom. However, as the reconstruction of the head phantom 
with noise shows, the algorithm according to equation 1 introduces artifacts along lines 
of high gradients of absorption values. This is well known from the expectation 
maximum likelihood (EM-ML) algorithm. 

The artifacts in the reconstructed image slice 31 occur as streaks. As may 
25 be seen from image slice 31 in Fig. 2, without any correction of the artifacts, i.e. 

without any artifact correction according to an exemplary embodiment of the present 
invention, a strong streak pattern can be observed. 

On the other hand, by reconstructing the image slice according to an 
exemplary embodiment of the method according to the present invention, in which the 
30 subtractions in the updates are weighted with the intrinsic statistical error <j Y} of the 
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9 

measured photon counts Y i during reconstruction of the image slice, the artifacts are 
compensated for and the streaks disappear: 

where a Yi = is the statistical uncertainty of a Poisson distribution. 

5 Each of the updates is hereby the result of a respective cycle of the iterative algorithm. 

It should be noted, that the method according to an exemplary 
embodiment of the present invention also works for a reconstruction with sub-sets of all 
measured projections. 

The image slices shown in Fig. 2 are the result of a convex algorithm 
10 (image slice 31), which is a special form of a ML algorithm, and a statistical uncertainty 
weighted improved convex algorithm (image slice 32), both after 500 iterations and 
starting from a uniform image. The level is 35 Houndsfield Units (HU), the window is 
100 HU. 

The convex maximum likelihood algorithm used for reconstructing the 
15 image slices of Fig. 2 is described in K. Lange and J. A. Fessler, "Globally convergent 
algorithms for maximum a posteriori transmission tomography", IEEE Trans. Med. 
Imaging 4, 1430-1450 (1995), which is hereby incorporated by reference. 

Fig. 3 shows a flow-chart of an exemplary embodiment of a method of 
artifact correction according to the present invention. After the start in step SI, an 
20 acquisition of the projection data set of the object of interest is performed in step S2. 
This may be performed by detecting a polychromatic x-ray beam emitted from a 
radiation source, which moves along a helical path around the object of interest. The x- 
ray beam passes the object of interest and is attenuated by object points inside the 
object of interest. During the movement of the radiation source a plurality of 
25 projections are acquired, each projection belonging to a respective mean position of the 
radiation source. After acquisition of the projection data set of the object of interest the 
number of detected photons Yi of an i-th projection is determined in step S3, e.g. with 
the help of a data processor. Then, in step S4, the statistical error a Y( of the number of 
detected photons Yi is calculated by the data processor and the vector p. of attenuation 
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parameters jlxj is set to an initial value, e.g. //,=!, in step S5. Each attenuation 
parameter (Xj belongs to a respective interval along the projection line of the i-th 
projection. 

In step S6 then, a counter is set to n = 0 after which the iterative loop for 
the image reconstruction starts. In step S7 the n-th update (which is in case of the 
beginning of the iterative loop the first update) of the attenuation parameter jlxj is 
calculated by 



T.Ucrl 



,,"+' - + — . 

' ' ' EA <>.,v 

Then, the method continues to step S8, where it is determined whether an 

10 end criterion is met or not. The end criterion may, for example, be that the iteration has 
been performed for a certain predetermined number of times or that a difference 
between consecutive updates does not exceed a certain predetermined threshold value. 
In case that it is determined in step S8 that the end criterion is not met, the counter is 
increased by one (step S9) and the method continues with step S7, where (n + l)-th 

15 update of the attenuation parameter fij is calculated. Step S7 - S9 may be iteratively 
repeated until the end criterion is met. In case it is determined in step S8, that the end 
criterion is met, the method continues to step S10, where it ends. 

It should be noted that, as may be apparent to the skilled person, the 
above described technique may be applied to all known reconstruction techniques in 

20 which a difference is summarized such as in the convex algorithm described above. 

Fig. 4 depicts an exemplary embodiment of a data processing device 
according to the present invention for executing an exemplary embodiment of a method 
in accordance with the present invention. The data processing device depicted in Fig. 4 
comprises a central processing unit (CPU) or image processor 151 connected to a 

25 memory 152 for storing an image depicting an object of interest, such as a patient. The 
data processor 151 may be connected to a plurality of input/output network or diagnosis 
devices, such as an MR device or a CT device. The data processor may furthermore be 
connected to a display device 154, for example a computer monitor, for displaying 
information or an image computed or adapted in the data processor 151. An operator or 

30 user may interact with the data processor 1 5 1 via a keyboard 155 and/or other output 
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devices, which are not depicted in Fig. 4. 

Furthermore, via the bus system 153, it is also possible to connect the 
image processing and control processor 1 5 1 to, for example a motion monitor, which 
monitors a motion of the object of interest. In case, for example, a lung of a patient is 
5 imaged, the motion sensor may be an exhalation sensor. In case the heart is imaged, the 
motion sensor may be an electrocardiogram (ECG). 



